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by 
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SUMMARY 

In  principle,  a  turbulence  field  is  to  be  governed  by  the  Navier-Stokes  equa¬ 
tions.  In  orde-  to  avoid  the  difficulty  of  treatment  due  to  the  non-linear  character¬ 
istics  of  the  Navier-Stokes  equations,  we  begin  with  the  assumption  that  a  turbulence 
field  may  be  represented  by  a  proper  distribution  of  many  elementary  vortex  lines, 
each  of  which,  being  a  particular  solution  of  the  Navier-StoKes  equations,  exhibits 
full  characteristics  of  the  non-linear  equation*.  Based  on  this  assumption,  we 
introduce  an  equation  which  governs  the  distribution  of  those  elementary  vortex  lines, 
in  the  same  way  as  the  Liouville  equation  governs  the  distribution  of  particles.  With 
respect  to  a  two-dimensional  field,  it  is  shown  that  Taylor' s  parabolic  correlation 
mode  for  short  distances  and  Kolmogoroff'  a  2/3-power  correlation  mode  for  moder¬ 
ate  distances  are  unified  to  one  correlation  moae  which  is  valid  for  the  entire  range 
of  correlation  distances.  With  respect  to  three-dimensional  case*,  introducing  re¬ 
marks  are  given  in  the  appendix. 
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I.  INTRODUCTION 


It  has  been  widely  felt  that  a  satisfactory  theory  of  turbulence  must  be 
one  in  which  non-linear  characteristics  of  fluid  motion  are  fully  taken  into  considera¬ 
tion.  It  is  the  purpose  of  this  paper  to  propose  a  theory  attempted  from  the  above 
viewpoint.  The  medium  is  assumed  to  be  incompressible,  and  detailed  treatments 
are  made  only  of  a  two-dimensional  case. 

As  is  well-known,  the  significance  of  velocity  correlation  as  the  cause  of 
transport  phenomena  in  a  continuum  was  pointed  out  by  O.  Reynolds  *  in  1895,  and 
clarified  by  G.  I.  Taylor  ^  in  19^1.  As  we  see  in  works  published  successively  by 
Taylor  Prandtl  4  von  Karman  Kolmogoroff  ^  and  others  between  then  and 
World  War  II,  earlier  treatments  were  mostly  kinematical  and  phenomenological, 
being  supported  by  similarity  principle;  physical  characteristics  of  turbulence  were 
described  in  terms  of  such  quantities  as  mixing  length,  size  of  eddy,  and  so  forth 
in  those  works.  The  nature  of  those  characteristic  lengths,  as  physical  quantities, 
was  known  only  vaguely  and  intuitively.  Among  those  achievements,  Taylor's 
experimental  results  of  velocity  correlation,  his  approach  of  coarse-graining  the 
Navier -Stokes  equations,  Kolmogoroff  s  correlation  functions  have  given  strong 
influences  on  modern  turbulence  studies. 

After  the  war,  Heisenberg,  Onsager,  Burgers,  Chandrasekhar  and 
7-8 

ethers  advanced  various  theories  in  attempting  to  find  feasible  statistical  laws 
governing  the  energy-spectrum  characteristics  of  turbulence.  The  approaches  of 
these  authors  were  inductive;  they  conceived  mutual  interactions  among  vortexes 
as  the  cause  of  thermalization  or  decay  or  irreversible  transfer  of  energy  among 
different  spectrums,  and  attempted  to  formulate  the  transfer  process  by  invoking 
statistical  methods  in  the  theory  of  Brownian  motion. 
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In  those  treatments  before  and  after  the  war,  we  see  that  those  authors 
were  assuming,  if  implicitly  and  intuitively,  that  turbulence  is  a  manifest  of 
vortexes  of  complex  distributions. 

Later,  particularly  in  the  last  ten  years,  attempts  have  been  made  by 
many  authors  (for  example, see  Refs.  9  and  10)  to  derive  feasible  laws  governing 
the  irreversible  process  in  a  turbulence  field  by  applying  statistical  treatments 
to  the  Navier -Stokes  equations.  Although  the  methods  vary  from  one  author  to 
another,  the  principle  is  one:  A  turbulence  field  may  be  described  precisely  in 
terms  of  the  velocity  correlation  factors  of  all  the  orders  (of  an  infinite  number). 

As  was  suggested  by  Taylor,  von  Karman  and  others  previously,  if  one  derives 
equations  governing  velocity  correlation  factors  from  the  Navier -Stokes  equations 
and  solves  them  properly,  one  may  achieve  his  goal.  In  this  venture,  however, 
one  must  meet  great  difficulties  which  are  common  in  dealing  with  non-linear 
phenomena.  Most  usual  methods  of  expansion  of  variables  in  series  are  not  feas¬ 
ible,  because  of  the  difficulty  of  convergency. 

We  have  had  many  experiences  of  overcoming  the  common  difficulty  of  non¬ 
linear  problems.  As  Eyring  **  wrote  recently,  science,  in  its  various  fields,  has 
invented  various  models  to  make  the  best  compromise  between  the  infinite  detail 
of  reality  and  the  limit  of  tractability.  It  should  be  remarked  that,  in  general,  a 
model  of  a  system  is  made  of  integrals  (invariants)  of  motion  of  the  system.  (  Of 
course,  approximate  integrals  of  motion  may  be  useful, too.  )  The  present  attempt  of 
turbulence  theory  is  made  also  in  this  sense. 

The  gists  of  the  present  attempt  are  1)  to  take  into  consideration  the  effect  of 
non-linear  characteristics  of  dynamical  processes  governed  by  the  Navier-Stokes 
equations,  and  at  the  same  time  £)  to  avoid  mathematical  difficulties  which  are  usual 


in  treating  significantly  non-linear  fields. 

In  other  words,  this  is  an  attempt  to  Lreat  the  structure  of  tu*Dulence  con¬ 
ceived  b>  Heisei  erg,  Onsager  and  others  in  a  deductive  sense.  The  consideration 
of  non-linear  characteristics  is  believed  to  be  essential  in  a  deductive  approach  to 
turbulence  theory,  particularly  when  the  turbulence  is  strong.  While  the  proposed 
deductive  approach  to  turbulence  may  appear  to  differ  from  those  of  other  workers, 
it  by  no  means  contradicts  or  dismisses  the  existing  understanding  of  turbulence  as 
accumulated  over  the  past  fifty  years.  Rather,  the  present  attempt  is  seen  to 
constitute  a  natural  and  rational  synthesis  made  of  th-.  Knowledges  achieved  by  the 
pioneers. 


II.  ELEMENTARY  VORTEX  LINES 


There  are  two  effects  which  characterize  a  flow  field  governed  by  the  Navier- 
Stokes  equations:  1)  Non-linearity  and  i)  stress  due  to  viscosity.  We  begin  our 
investigation  of  turbulence  by  ignoring  the  effect  of  viscosity.  This  approximation 
is  valid  if  the  characteristic  Reynolds  number  of  turbulence  is  sufficiently  large, 
insofar  as  turbulence  is  investigated  in  a  free  flow.  (See  section  VI)  On  this  ap¬ 
proximation,  we  have  to  make  an  important  choice  between  the  following  two  views: 

1.  "A  turbulence  field  is  a  field  of  vorticity,  continuously  but  non-unif  rmly 
distributed  in  the  configuration  space.  The  strength  of  vorticity  is  finite.  There¬ 
fore,  there  is  no  discrete  vortex  line,  or  tube  conceived  by  Helmholtz.  Based  on 
this  consideration,  a  turbulence  field  is  to  be  treated  as  a  continuous  field  of 
vorticity  governed  by  the  Navier-Stokes  equation". 

i.  "Due  to  initial  and  boundary  conditions,  such  as  an  array  of  rods  inserted 
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in  a  flow  which  is  otherwise  uniform  ,  however,  aggregations  of  vorticity  may 
exist,  ihe  field  of  vorticity  may  be  continuous  in  the  precise  sense.  But  in  an 
approximation,  it  is  possible  to  represent  the  field  as  an  ensemble  of  discrete 
vortex  lines  or  tubes.  By  considering  kinetic  theory  of  such  discrete  vortexes,  it 
may  be  possible  to  obtain  essential  characteristics  of  a  turbulence  flow,  even  if 
we  may  be  ignoring  the  microscopic  and  detailed  structure  of  the  flew  field.  ” 

As  supporting  the  second  approach,  we  have  the  following  two  knowledges: 

A  vortex  line,  which  may  approximate  a  local  aggregation  of  vorticity,  is  a  particular 
solution  of  the  non-linear  Navier-Stokes  equations,  or  the  Euler  equations  on  ignoring 
the  viscosity  effect,  and  is  invariant  as  governed  by  Helmholtz1  s  three  laws. 

Secondly,  we  have  a  successful  experience  of  kinetic  theory  of  gases  on  the  assump¬ 
tion  that  a  gas  consists  of  discrete  and  stable  particles  in  spite  of  the  fact  that  the 
most  precise  description  of  a  gas  as  a  whole  is  a  quantum -mechanical  wave  function, 
which  is  spread  in  the  entire  space.  Treating  a  gas  as  described  by  a  wave  function 
might  appear  precise,  but  in  fact  we  may  miss  its  essential  characteristics  due  to 
its  enormous  complexity.  Instead,  the  approximate  model  of  particle  structure 
enal'es  us  to  see  main  characteristics  of  a  gas  rather  easily.  This  is  the  approach 
of  kinetic  theory. 

According  to  the  above  consideration,  we  assume  that  a  turbulence  field  is 
represented  by  an  appropriate  distribution  of  elementary  vortex  lines;  each 
elementary  vortex  line  is  of  the  same  intensity  and  behaves  according  to  Helmholtz' s 
three  laws.  A  turbulence  field  consists  of  such  elementary  vortex  lines  distributed 
non-uniformly.  A  vortex  tube  of  a  large  size  may  be  composed  as  a  cluster  of  many 
elementary  vortex  lines  of  the  same  direction.  It  is  noted  that  such  a  vortex  tube 
of  a  large  size  exhibits  a  sort  of  thermalization  phenomenon  due  to  random  migrations 
of  elementary  vortex  lines  constituting  the  tube,  as  is  illustrated  in  Fig.  1. 
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Fig,  - 1  -  von  Kirman'  a 


vortex  street: 

Schematic  representation 
of  the  vortex  street  by 
means  of  elementary 
vortex  lines. 
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a  turbulence  flow  behind 
an  array  of  rods  inserted 
in  a  uniform  flow. 

Each  eddy  consists  of 
many  elementary  vortex 
lines. 
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The  distribution  of 
elementary  vortex  lines 
near  the  final  stage  of 
thermalization  of  the 
turbulence  produced  in 
(b) 
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ill.  KINETIC  EQUATIONS  OF  ELEMENTARY  VORTEX  LINES 

We  aasume  that  a  turbu'ence  field  consists  of  N  elementary  vortex  lines  which 
are  parallel  to  th*  ?. -axis  and  of  which  one  half  is  of  the  positive  sign  and  the  other 
of  the  negative  sign.  We  may  define  the  distribution  of  the  N  vortex  lines  by 
D(N)=n6(r,  -r*)6(5.  0*)  (3.1) 


where  r.  is  a  function  of  time  and  denotes  the  position  of  vortex  i,  and  fl.  denotes  ti 

$ 

direction  and  intensity  of  vortex  i.  Of  course.  f).  is  an  invariant  in  the  present  case, 
r.  and  Oj  are  independent  variables.  If  vortex  CL  is  at  r^  and  induces  velocity  v. ^ 


at  r., 


-  _  (ri-rj)  x  nj 

i‘i  (r.-r.)^ 

'  i  j' 


(3.<i) 


(3.4) 


—  ‘  J 

If  V.  is  the  velocity  of  the  main  (average)  flow, 

t  -•  _  -a 

vvi + 1  vij  <»•  ■» 

is  the  total  velocity  at  rj :  the  trajectory  of  vortex  i  at  is  governed  by 

dri  *  V 

W  ‘  (3.4) 

Therefore  D^def*ned  by  (3.  I‘,  atisfies 

(|t+y  V  4  JD^O  (3.5) 

ot  ^  1 

i  i 

The  above  equation  is  similar  to  the  Liouville  equation  of  N  particles. 

In  a  manner  similar  to  that  of  the  ordinary  kinetic  theory  of  particles,  we  may  define 

FV,)=  oW-.  to;  d', 


FW(‘-  J)  *  D(NkJi.  j  drA 

and  so  forth.  Consideration  of 

L. a  x^i 

%?i'Vij  "J?i  '  (r>V 


leads  to 


r ...  .  X-  (Ll<[v.  =  o 

J  ^rj  J*,5|rj  lJ  3 


(3.  7) 
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(3.8) 


Hencj,  it  is  easily  shown  that  Eq.  (3.  5)  is  reduced  ‘o 
(  2_+  v  3  .  i  ,  \  •  T.  *  . 


Kt+Vi  >  F  F  (ijJd-Y 


where 


dx.=dr.  do. 
J  1  J 


[|r+  +V  •  >?-*  *  V%  ] F  w«> 

*1  (vik  •  i^r  > F<J,(i>  j- 10  d  V° 

w  1 


j  [=1,  i.,  ....  N 


(3.  10) 


and  so  forth.  These  equations  constitute  an  indefinite  series  of  equations  tn  a 
manner  similar  to  that  of  the  Bogoliubov-Born-Green  Kirkwooc  Yvon  hierarchy 
of  equations  governing  the  distributions  of  the  molecules  c  ostituting  a  gas. 


IV.  MODES  OF  INTERACTION  AMONG  ELEMENTARY  VORTEX  LINES 
1.  Microscopic  Interaction.  If  the  interaction  between  elementary  vortex 
lines  i  and  j  is  strong,  and  they  are  fairly  remote  from  the  other  vortexes,  the 
interaction  is  almost  binary.  In  this  case,  Eq.  (3.10)  yields 

$+V,-*-+Vst: 

»  J  _  - 

(V.=V.-0  is  assumed)  (4.  1) 

1  j 

According  to  the  above  equation,  two  vortexes  i  and  j  rr. ove  along  trajectories 


such  as  illustrated  in  Fig.  i. 


/  ~  \ 
«  t 

V  C. 


Fig.  i  -  A  pair  of  elementary  vortex  lines 
rotates  if  their  signs  are  the  same, 
but  they  move  on  two  parallel  lines 
if  their  signs  are  different. 


’J 


If  i  and  j  are  of  the  same  sign,  they  move  on  a  circle,  while,  if  their  signs 
are  mutually  opposite,  they  move  on  two  parallsl  lines.  Their  interaction  does  not 
terminate  unless  they  are  disturbed  by  other  elementary  vortex  lint  .  We  assume 
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that  the  strong  correlation  between  i  and  j  terminates  in  a  finite  and  microscopic 
time  period  due  to  weak  but  frequent  pertubations  due  to  ve»  cities  exerted  by  other 
vortexes  such  as  represented  by  the  last  integral  term  in  Eq.  (J.  11).  According  to 


(4. 1),  we  obtain 


/re  obtain 

,-i- F«,i, j)  d  x.  -j .  -i- 1 . 4_  ,  j,  d,. 


-f<4r  ‘  j)dx; 


J 


(4.  a) 


where  d/dt  is  an  operator  meaning  the  differentiation  along  the  trajectories  of  the 
two  vortexes  i  and  j.  Since  it  is  shown  that 


is  negligible,  we  obtain 


W 

(i.  j)  dx  )  dt 
j 


.—  F^)(i,  j)  dx\  dt 

ij  j 

=  T~JJi t  F<^(i«j>dxj  dt 

=  J  [  F'^tM-F^oJdx 

=  J  [  F(1)(r.  +  :r)  T'Xhr.^r.) 

-F(n(r.)  F(l,<r.)  ]dx  (4.J) 

Here  i  r.  and  1  r .  are  the  distances  traveled  respectively  by  vortexes  i  and  j  during 
a  continuous  correlation  period  of  time  r.  The  approach  deriving  Eq.  (4.  j)  and  the 
interpretation  are  the  same  as  those  employed  for  the  Boltzmann  collision  integral 
derived  from  the  BBGKY  hierarchy  under  the  binary  collision  assumption,  except 
that  1)  the  integrand  in  Eq.  (4.  J)  is  due  to  the  non-uniformity  of  the  vortex  distribu¬ 
tion  in  the  configuration  space  while  the  integrand  of  the  Boltzmann  collision  integral 
is  due  to  the  non-uniformity  of  the  particle  distribution  in  the  momentum  space,  and 
that  i)  the  termination  of  the  present  interaction  between  vortex  ;  and  vortex  j  is 
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caused  by  disturbances  of  other  vortexes,  while  the  termination  of  a  collision  between 
two  particles  is  due  to  the  nature  of  Newtor's  dynamics  and  the  relevant  force  law. 

(In  case  of  charged  particles,  the  termination  of  a  binary  collision  is  due  to  dis¬ 
turbances  of  other  particles  ) 

L,  interactions  Among  Semi-Stable  Clusters  of  Elementary  Vortexes. 

A  two-dimensional  turbulence  field  is  produced  by  an  array  of  rods  placed  perpen¬ 
dicularly  in  a  flow  which  is  otherwise  uniform.  In  origin,  each  rod  produces  a 
street  of  vortexes  similar  to  von  Harman*  s  vortex  street.  Of  course,  there  are 


interactions  among  those  streets  formed  by  many  neighboring  rods.  Each  vortex 
constituting  the  vortex  streets  is  considered  as  a  cluster  of  elementary  vortex  lines 
and  may  be  initially  simulated  by  a  vortex  tube  of  finite  diameter  instead  of  a  vortex 
line,  (dee  Fig.  1).  Let  us  suppose  that  two  vortex  tubes  are  in  interaction.  Because 
of  the  finiteness  of  the  diameter  of  a  tube,  their  interaction  is  more  complex  than  an 


interaction  between  two  vortex  lines 

Fig.  3 


A  is  a  cluster  of  elementary  vortex  lines  of 
negative  sign,  and  B  is  of  positive  sign.  Initially, 
their  cross-sections  were  circular.  As  they  move 
on  two  parallel  lines,  their  cross-sections  are 
deformed:  The  part  of  A  closer  to  B  proceeds  faster 
than  the  other  part  of  A,  etc. 


In  Fig.  3,  two  vortex  tubes  of  different  signs,  A  and  B,  are  in  interaction.  They  are 
mutually  driven  on  two  parallel  lines  in  the  first  approximation,  but  the  part  of  vor¬ 
tex  A  remote  from  B  proceeds  slower  than  the  part  close  to  B.  Similarly,  the  part 
of  B  remote  from  A  proceeds  slov.er  than  the  part  close  to  A.  Furthermore,  each 
vortex  tube  is  rotating  by  itself.  As  a  result,  the  cross-section  of  each  vortex  tube 
is  deformed  as  time  passes.  Similar  and  more  complex  deformations  are  caused  on 
a  vortex  tube  surrounded  by  many  ot  er  vortex  tubes.  Thus  a  cluster  (tube)  decays 
from  its  outer  fringe  toward  its  center.  As  a  cluster  decays,  those  elementary 
vortexes  which  have  drifted  outward  may  be  mixed  with  elementally  vortexes  of  the 
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opposite  sign  which  have  drifted  from  another  cluster.  Those  two  groups  of  elemen¬ 
tary  vortexes  of  different  signs  may  be  mixed  at  random,  and  in  effect  the  diameter 
of  a  vortex  cluster  becomes  smaller.  As  this  process  of  decay  proceeds,  each 
cluster  becomes  thinner  and  the  effect  of  each  cluster  deforming  other  clusters  be¬ 
comes  weaker.  This  situation  is  conceivable  if  we  notice  that  a  vortex  line,  instead 
of  tube,  is  not  deformed  at  all  by  other  vortex  lines.  This  may  be  the  reason  for  the 
existence  of  an  almost  homogeneous  and  steady  turbulence  field  indifferent  of  its 
initial  method  of  formation.  According  to  the  above  consideration,  however,  there  is 
no  homogeneous  and  steady  turbulence  in  the  strict  sense.  A  seemingly  homogeneous 
turbulence  is  still  on  the  process  of  decay,  although  the  rate  of  decay  may  be  small. 

3.  Interactions  of  the  Vlasov  Type.  If  a  distribution  of  elementary  vortex 
lines  is  non-uniform  in  the  microscopic  sense,  the  integral  term  in  Eq.  (3.  8)  does 
not  vanish  even  when  the  microscopic  correlation  between  i  and  j  is  ignored;  that 

y  V  F(1)(j)dx  .4  F(1)<i) 

'  'J  J  on 

j 

does  not  vanish  and  gives  the  effect  of  macroscopic  interaction. 


V.  CORRELATION  OF  VELOCITY  FLUCTUATIONS 
DUE  TO  CLUSTERS 

As  is  described  in  the  last  section,  clusters  of  elementary  vortex  lines  are 
assumed  to  constitute  a  homogenous  turbulence  field.  Based  on  this  model,  it  would 
be  possible  to  calculate  the  correlation  between  fluctuations  at  two  different  positions 
in  the  space.  Such  correlations  were  first  observed  experimentary  by  G.  I.  Taylor 
in  1935.  A  schematic  correlation  curve  is  given  in  Fig.  4. 
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Fig.  4-1;  Parabolic  law  by  Taylor;  II;  III 
power  law  by  Kolmogoroff;  the 
present  theory  gives  the  solid 
line. 


*  fi 


10 


Taylor  could  derive  the  correlation  curve  from  the  Navier-Stokes  equations  for 
distances  between  two  positions  which  are  short  as  compared  with  the  size  of  eddy; 
that  is  the  parabolic  law: 

vAx  vBx-o(l-R^|  j 

Later,  Kolmogoroff  gave  a  formula  of  correlation  valid  for  a  broader  range  of  dis¬ 
tances  of  two  positions  where  the  correlation  is  considered,  by  means  of  a  dimension¬ 
al  consideration;  that  is 

Rils'  11 

but  so  far,  no  rational  attempt  has  been  successful  in  deriving  the  correlation 
function  which  is  valid  over  the  entire  domain  of  R.  In  the  following,  it  will  be  show 
that  our  present  model  is  useful  lor  the  derivation. 

We  assume  that  each  cluster  consists  of  n  elementary  vortexes  and  is  a 
circular  cylinder  of  diameter  d;  one-half  of  the  number  of  the  clusters  is  of  the  posi¬ 
tive  sign  and  the  other  half  of  the  negative  sign.  They  are  distributed  at  random 
in  the  configuration  space  with  a  uniform  density.  We  also  assume  that  the  measure¬ 
ment  of  velocities  is  made  with  respect  to  a  coordinate  axis  system  which  is  moving 
with  the  ''elocity  of  the  main  (average)  flow.  First.it  will  be  shown  that  the  velocity 
fluctuation  at  a  position  is  attributed  mainly  to  the  nearest  neighboring  vortex  cluster. 
Secondly,  we  calculate  the  velocity  correlation  at  a  pair  of  positions. 

1.  Localization  of  the  Random  Pstribution  of  a  Cluster.  As  discussed  pre¬ 
viously,  a  cluster  of  elementary  vortex  1  nes  is  produced  by  a  proper  boundary  and/or 
initial  conditions  .  Since  the  velocity  of  a  cluster  is  finite  and  the  time  period  in 
which  a  cluster  passes  through  the  field  domain  under  investigation  is  also  finite,  the 
probability  distribution  of  the  cluster  is  localized.  The  local  domains  of  space  in 
each  of  which  a  cluster  is  present  with  uniform  probability  density  are  patched  to¬ 
gether  and  constitute  the  entire  space  domain  under  consideration.  This  consideration 
may  be  feasible,  if  the  formation  of  the  turbulence  field  is  done  in  a  proper  way;as  such 
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it  le  made  with  an  array  of  rods.  The  linear  dimension  of  such  a  local  domain  is 
denoted  with  which  is  of  the  order  of  the  distance  of  two  neighboring  rods.  By 
denoting  the  number  density  of  clusters  with  n,  n”*^is  of  the  order  of  the  distance 
between  two  neighboring  clusters.  Assuming  that  an  observer  is  moving  with  the 
velocity  of  the  main  flow  V,  and  considering  that  the  fluctuating  part  of  velocity 
of  a  cluster  v  is  smaller  than  V,  we  have 


L.  (the  distance  which  the  observer  travels  during 
his  observation) 


Hence 


l- 


l 

v 


L  >  >  n 


■Ml 


Lv 


v  -Ml 
>  >  —  n 

V 


Since 


—  <  1  ,  we  may  assume  that 
V  -1/2 


(5.  1) 


is  the  condition  of  our  observation.  This  relation  may  be  given  based  on  von  Karman's 
vortex  itreet  theory.  This  is  also  the  condition  that  an  arbitrary  position  in  the  field 
has  the  same  probability,  as  the  averaged  one,  to  be  occupied  by  a  cluster.  A  similar 
condition  is  realized  in  the  special  case  where  each  cluster  is  distributed  with  a  uni¬ 
form  probability  density  all  over  the  field.  In  both  cases,  the  probability  densities 
of  clusters  are  the  same.  With  respect  to  fluctuations  in  density,  however,  the  two 
cases  are  different:  In  the  former  case,  a  space  domain  of  linear  dimension  larger 
than  t  cannot  be  completely  void  of  any  one  cluster  at  a  moment  of  time.  On  the 
other  hand,  in  the  latter  case,  it  is  even  possible  that  the  entire  field,  except  for 
one  spot  at  least,  is  void  of  any  cluster. 

1.  Significance  of  the  Nearest  Neighboring  Cluster. 

The  most  significant  contribution  of  velocity  fluctuation  at  a  point  p  is  made  by 
the  cluster  which  is  the  nearest  neighbor  of  p.  The  reasons  are: 


ll 


1.  The  nearest  distance  and  hence  the  largest  induced  velocity. 
i.  The  location  of  a  cluster  is  localized  even  though  the  precise  location 
at  each  moment  is  unknown.  We  should  notice  the  difference  between  the  uniform 
probability  distribution  of  a  cluster  and  the  localized  probability  distribution  of  the 
cluster.  To  explain  the  difference,  we  suppose  that  cluster  a  of  strength  w  is  local¬ 
ized  in  domain  A  and  the  cluster  b  is  in  domain  B.  If  they  are  of  the  same  sign  and 
of  the  same  strength,  and  A  and  B  are  mutually  symmetric  with  respect  to  P  where 
induced  velocities  are  observed,  the  order  of  the  fluctuation  induced  by  the  two 
clusters  at  P  is 

JS _  8 
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where  L  9  is  the  angle  with  which  A  and  B  is  seen  from  P,  and  L  is  the  approximate 
distance  of  A  and/or  B  from  P.  On  the  other  hand,  if  they  are  completely  at  random 
with  respect  to  direction,  the  order  of  fluctuation  is  to  be  of  the  order  of 

in 

irr'L  (5. 3) 


In  order  to  investigate  the  order  of  the  magnitude  of  fluctuations  at  P  induced 


by  clusters  surrounding  the  point,  we  assume  that 


rr  ■£, 


L  . 


is  of  the  order  of  (density 


of  clusters)"*.  Also,we  assume  that  tt  is  the  area  of  the  domain  in  which  a  cluster 


moves  at  random  during  the  time  period  of  our  investigation. 

Fig.  5  -  The  space  is  divided  with  co-centric  circles  of 
~  which  the  center  is  at  P,  where  the  velocity 
fluctuation  is  observed,  and  the  radii  are 
3<s  5 1,  ... 


We  divide  the  space  with  co-centric  circles  of  which  the  center  is  at  P  and  the 
radii  are  it,  5{.,  . . .  The  domain  of  the  smallest  circle  is  named  domain  1. 

Domain  i  is  the  domain  encircled  by  the  smallest  circle  and  its  neighboring  circle. 
Similarly  we  may  define  outwardly  domain  3,  domain  4,  . , .  ,  domain  5,  . .  . 

Domain  i, which  is  between  the  circle  of  radius  {  and  radius  3/  is  divided  by  radial 
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lines  into  cells  of  the  same  area  rr  •l/'  which  is  the  area  of  domain  1.  The  number  of 
the  cells  is 


- I - 

n-r 


-  8. 


Suppose  that  a  cluster  of  the  same  sign  is  confined  in  each  cell.  The  velocity  fluc¬ 
tuation  at  P  induced  by  the  cluster  in  domain  1  is  of  the  order  of 


V  _  W 

1  "  7ttT7 T- 


(5.4) 


The  total  fluctuation  of  velocity  at  P  induced  by  the  eight  clusters  in  domain  i  is  of 
the  order  of 


where 


*  IrtT -  A  03  A  e  3 )  <4  L  0  a  ) 

3 

L,  =1  l,  L  0a  =  =  -J-  . 


Va  =0.  11  X  Vy 


It  is  easily  shown  that 

(5.5) 

Similarly  we  may  consider  the  effect  of  the  clusters  in  domain  3.  The  number  of  the 
cells  is 


n(5 4J3-'1  (5lP  _ 


l  a 


=  16 


and 


where 


Hence 


V3= 


w 


i  ttA 


(^A93)  (4A03)  (8A03) 
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La=  ’  A0a=  mi  xT5  “"TB' 


1,1.4  v- 

v3=  v3  7(7}  x  8  =  T‘  (5.6) 

In  general,  the  effect  of  the  clusters  of  the  vth  domain  is  given  as  follows:  The  number 


of  the  cells  is 


and 


-ill  lilgi!  .  8  (v.„ 

TT  l* 


49  • 


v  -lx  8(V-1)  8  (V-1)  , 


*  (v-D  4 


Hence 


v  a. 

v 


w 


<dnLv 


(  A  9v)  \>  1,  i, . . . ,  i  v 


-1 


(5.7) 
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Here 


(5.8) 


<2kv=8  (v  1) 

as  a  function  of  v  is  tabulated  as  follows  ; 

v  <23456789 
kv  3  4  4<5  5  5^6  5<6'  5  6"  6 

where  5<6,  for  example,  means  a  value  between  5  and  6.  It  is  an  easy  matter  to 
see  that 

vi  va  ’  va  ■ 

constitute  a  series  which  converges  very  quickly,  and  v^  is  a  good  approximation 
of  the  entire  fluctuation. 


3.  Correlations  of  Velocities. 

In  order  to  see  consequences  of  the  above  consideration,  we  calculate  cor¬ 
relations  between  velocities  at  two  positions  A  and  B,  induced  by  cluster  C  which 
is  the  common  nearest  neighboring  one  of  the  two  positions  A  and  B.  It  is  assumed 


that  A  and  B  are  moving  with  the  main  flow. 


Fig.  6  -  We  calculate 

correlations  between 
velocities  at  two 
positions  A  and  B, 
induced  by  Cluster  G 


The  distance  between  A  and  B  is  R.  R  is  assumed  to  be  changed  from  zero  to  in¬ 
finity.  If  R  is  larger  than  i,  the  nearest  neighboring  cluster  of  A  is  different  from 
that  of  B.  Then,  in  the  approximation  of  ignoring  those  clusters  other  than  the 
nearest  neighboring  one,  there  is  no  correlation  between  the  fluctuations  at  A  and  B. 
We  also  assume  that  the  radius  o  of  a  cluster  is  finite  and  is  sufficiently  smaller 
than  l.  A  and  B  are  on  the  x-axis  with  the  origin  of  the  coordinate  system  at  the 
center  of  AB.  It  is  often  convenient  to  give  the  position  of  C  in  terms  of  r,  the 
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distance  from  the  origin  O  and  0,  the  direction  angle  with  respect  to  the  x-axis. 
Also  we  define 


pA=str. 


0  =<:CAO, 
A 


Pb:BC- 


g  =n.^CBO. 


1.  The  velocities  at  A  and  B. 


(i)  If  s^*'b  >  O*  we  have  for  the  velocities  at  A  and  B  induced  by  cluster  C 
of  intensity  uj 
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<*ttP 


A 

cos  9 a 
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cos  0 
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It  is  a  simple  matter  of  geometry  to  find 


=  [  A  Rr  cose]  ^ 1 

B=f  (-T—  )^+  r'5-  R  r  cos  o] 

1  '  R  J 


sin  9  = 


r  sin  9 


A  n 


sin9B= 


A 

r  sin® 

°B 


cos  9  A  = 
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\U 

+  r  cos  0 


cos  eB= 


+  r  cos  0 
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(5.9) 
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It  is  r 


ly  seen,  in  view  of  the  symmetry  of  the  field,  that 
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Hence 
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3.  v,  v. 
Ax  Ay 


It  is  easily  shown  that  this  correlation  factor  vanisl  es: 


v  v* 
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=  O. 


4.  v.  /_ 
Ax  Bx 


In  order  to  calculate  this  factor,  vj^xvgx  to  be  integrated  over  the 
domain  of  r,  in  which  a  point  is  within  distances  less  than  /  both  from  A 
and  B.  The  manipulation  may  be  complicated,  depending  on  the  relations 
among  R,  c  and 

(i).  If  f-  <R,  A  and  B  do  not  have  any  common  nearest  neighboring 
clus^.r.  Hence 


v 


AxvBx 


=  O 


(5.  15) 
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The 


(ii).  If  R  >  i,  R  -  6  >  CT,  it  is  easy  to  calculate  vi^xvgxi 
integration  is  to  be  made  over  the  domain  abed  indicated  in  Fig.  The 
maximum  value  of  1  for  given  6  is  the  solution  of 


or 


r3  +  (4-)  +  i  (^r~)  r  cos  9=^ 
m  c  c  m 

3  8-  1)1 

)  o  i  •  J 


rm~ 


8  6 


(5.  16) 


Vre  take  only  the  positive  value  of  r 


m 


Fig.  7-  If  R  >  l,  the  clusters  only  in  the  domain 

abed  contribute  to  the  correlation  between 
A  and  B. 


Then  we  have 


v .  v_ 
Ax  Bx 
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r  sin  or  d  r  d  8 
[-jp-  +  r*J3  -R^-  acos: 
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Since  R  is  always  larger  than  r,  we  expand  the  integrand  in  powers  or  r/R. 
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Noting  that 
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fr  sin2  9  (  -  JL  +  16  —  ~.9 ....  )  drd9 
I  R3  R3 
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we  obtain 
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(iii)  R  <<  o  <1 •  The  condition  is  illustrated  in  Fig.  8.  A  and  B 
have  a  common  nearest  neighboring  cluster.  The  domain  where  r  <  a  is 


called  domain  I. 


Fig.  8  -\iP«  ^ ,  A  and  B  have 
a  common  nearest 
neighboring  cluster. 


If  C  is  in  domain  I,  we  have 
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If  C  is  inside  the  domain  II  where  r  <o.  A  and  B  are  present  within  the 
core  of  a  cluster.  Hence, 


(VAxVBx )= 


^^rdrde 
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Therefore,  we  readily  have 
,  a 
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Fig.  -  4  -  I;  Parabolic  law  by  Taylor; 
II;  i.1  3  power  law  by 
Kolmogoroff;  the  present 
theory  gives  the  solid  line. 


VAxvBx  8^ven  by  (5.  18)  is  valid  for  R  <0  and  is  in  agreement  with  the 
experimental  result  obtained  by  Taylor  (Fig.  4,  a  )t v^xvgx  given  by 
(5.  17)  is  valid  for  R  >  l  +  a,  and  ’8  in  agreement  with  the  part  y  in  Fig. 

4.  It  is  naturally  expected  that  part  0  in  the  same  figure  must  be  obtained 
if  we  carry  out  the  calculation,  with  no  particular  difficulty,  by  a  means 
of  a  compu'ng  machine. 
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VI.  DISCUSSIONS  AND  CONCLUDING  REMARKS 


1.  The  common  difficulty  in  ordinary  approaches  of  turbulence  theory 
is  *hac  equations  of  correlation  functions  derived  from  the  Navier-Stokes  equa¬ 
tions  t  re  not  closed.  The  difficulty  becomes  serious  due  to  the  fact  that  the 
Navier-Stokes  equations  are  non-linear.  In  order  to  have  a  closed  set  of  equations 
or  correlation  functions,  it  is  necessary  to  ignore  some  terms  in  some  equations. 
In  general,  however,  there  is  no  assurance  that  those  ignored  terms  are  really 
ignorable  without  any  significant  errors. 

By  noticing  that  a  vortex  line  is  a  singular  sc.ation  of  the  Navier-Stokes 
equations  where  viscosity  is  ignored  and  that  a  vortex  line  is  invariant  in  the 
sense  of  Helmholtz,  we  assume  that  a  turbulence  field  constitutes  many  similar 
(elementary)  vortex  lines.  On  this  assumption,  we  regard  a  turbulence  field  to 
be  a  discontinuous  field  of  velocity.  By  so  doing,  however,  the  non-linearicy 
of  the  field  is  tractable. 

3.  The  effect  of  viscosity  is  ignored.  One  should  note,  however,  that  the 

effect  of  viscosity  of  a  vortex  line  is  not  to  extinguish  any  amount  of  vorticity; 

the  vorticity  concentrated  on  an  elementary  vortex  line  diffuses  finally  due  to 

viscosity,  but  the  total  vorticity  is  invarJant.  If  this  broadening  of  a  vortex  line 

due  to  viscosity  does  not  seriously  affect  the  main  characteristics  of  the  field 

under  consideration,  it  would  be  feasible  to  ignore  the  viscosity  effect.  The 

situation  depends  on  the  time  and  space  scales  on  which  a  field  is  considered. 

It  is  well-known  that  the  distribution  of  vorticity  initially  concentrated  on  a  filament 

is  given  by  „  =JA_  exp  (  -jf.  ) 

t  4vt 


where  vis  Viscosity  coefficient)  /  (fluid  density),  r  the  distance  from  the  filament 
and  t  the  time.  (  See  L.  Prandtl,  the  Mechanics  of  Viscous  Fluids,  in  Aero¬ 
dynamic  Theory  edited  by  W.  F.  Durand,  Vol.  Ill,  p.  68.  )  According  to  this  for¬ 
mula,  the  order  of  the  broadening  of  a  filament  is  given  by 

r  =  (vt)  lU  (6.  1) 

If  r  given  here  is  sufficiently  smaller  than  the  linear  dimension  of  our  close 
observation  our  assumption  of  neglecting  viscosity  would  be  feasible.  In  this 

case,  it  would  be  proper  to  take  for  t  the  duration  of  our  observation;  that  is, 

l 

t  =  ,  t  >  r 

v 

where  v  is  the  intensity  of  velocity  fluctuation.  Substitution  of  the  above  in 

(6.  1)  leads  to  — 

'  f,v 

-  >  1.  (b.i) 

v 

Our  present  treatment  made  by  ignoring  viscosity  is  feasible,  if  (6. 1)  is  satisfied. 

4.  In  a  three-dimensional  turbulence  field,  the  behavior  of  a  vortex  line 
is  complex  and  the  equation  of  the  distribution  of  vortex  line  ■  s  also  complex  as 
is  shown  in  the  Appendix. 

5.  The  present  article  is  merely  a  conceptual  introduction.  The  treatment 
of  the  basic  equation  introduced  here  may  vary  significantly  depending  on  the 
characteristics  of  boundary  and  initial  conditions  with  which  a  turbulence  field 
exists.  The  situation  may  be  similar  to  that  of  the  Boltzmann  equation  in  cases 
of  rarefied  gas  dynamics. 

6.  We  note  that  there  is  a  resemblance  between  the  present  theory  and 
the  ordinary  kinetic  theory  of  particles.  Specifically,  Eq.  (3.  5)  is  similar  to  the 
Liouville  equation  of  many  particle  systems.  In  spite  of  the  similarity  of  forma¬ 
lism,  there  is  a  significant  difference  between  the  present  theory  and  the  ordinary 
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kinetic  theory  of  particles.  At  the  kinetic  stage  of  coarse-graining,  before 
fluid-mechanical  stage,  the  effect  of  interaction  among  vortexes  result  in 
"diffusion"  in  the  configuration  space.  On  the  other  hand,  the  effect  of  inter¬ 
actions  among  molecules  causes  friction  and  diffusion  in  the  momentum  space, 
as  is  well  seen  in  a  Fokker -Planck  type  equation  which  may  be  derived  from  the 
Boltzmann  equation. 

7.  It  should  be  noted  that  a  field  of  vorticity  is  not  always  a  field  of  turbu¬ 
lence;  it  is  possible  for  a  field  of  vorticity  to  be  stable.  Phenomena  in  a  stable 
field  of  vorticity,  for  example,  laminar  boundary  layer,  are  out  of  the  scope  of 
the  present  investigation.  Our  investigation  begins  with  the  assumption  that  a 
field  under  consideration  is  turbulent. 
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APPENDIX 

THREE-DIMENSIONAL  TURBULENCE 


In  a  three-dimensional  turbulence  flow,  an  elementary  vortex  line  is  not  a 

straight  line;  it  forms  often  a  closed  line  as  is  schematically  illustrated  in  Fig. 

9.  Let  us  cut  the  closed  line  into  segments  of  equal  length  l. 

Fig.  -  9  -  A  three-dimensional 
vortex  line  may  be 
closed.  We  cut  the 
closed  line  into  seg¬ 
ments  of  equal  length, 
and  investigate  the 
distribution  of  those 
elementary  segments. 

We  investigate  the  distribution  of  elementary  segments.  The  strength  of  an 

elementary  vorfjx  segment  is  invariant,  but  the  direction,  the  length  l  and  also 

the  curvature  change.  If  there  is  no  large  non-uniformity  of  veS^tLy  in  the  flow, 

l  may  be  almost  invariant;  if  the  segment  is  sufficiently  short,  we  may  ignore 

the  curvature.  By  assuming  these  we  have 


(  ri  “  rJ)  x 

-•  -♦a 


(  r.-r.) 

i  y 


As  similar  to  Eq.  (3.2)  here  r.  is  the  position  of  the  center  of  elementary  segment 
i  and  r^  of  elementary  segment  j;  *  is  a  function  of  i  which  accounts  for  the  finiteuess 
of  length  i.  Then  we  have  for  the  basic  Liouville  equation 


as  similar  to  Eq.  (3.  5). 
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The  last  term  in  the  left  hand  aide  ia 
elementary  segments;  note  that 

iLi_ 

d  t  l  i  1 

Aa  ia  illuatrated  in  Fig.  9>  elementary  aegment  i  and  its  neighboring  aegment  j 
which  belong  to  the  same  vortex  line  have  a  correlation.  This  situation  may  be 
accounted  for  effectively  by  the  consideration  that  the  probable  distribution  of 
a  segment  ia  localized. 
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13  ABSTRACT 

In  principle,  a  turbulence  field  is  to  be  governed  by  the  Navier- 
Stokes  equations.  In  order  to  avoid  the  difficulty  of  treatment  due  to  the 
non-linear  characteristics  of  the  Navier-Stokes  equations,  we  begin  with 
the  assumption  that  a  turbulence  field  may  be  represented  by  a  proper  dis¬ 
tribution  of  many  elementary  vortex  lines,  each  of  which,  being  a  particular 
solution  of  the  Navier-Stokes  equations,  exhibits  full  characteristics  of  the 
non-linear  equations.  Based  on  this  assumption,  we  introduce  an  equation 
which  governs  the  distribution  of  those  elementary  vortex  lines,  in  the  same 
way  as  the  Liouville  equation  governs  the  distribution  of  particles.  With 
respect  to  a  two-dimensional  field,  it  is  shewn  that  Taylor’ s  parabolic  cor¬ 
relation  mode  for  short  distances  and  Kolmogoroff’ s  2/3-power  correlation 
mode  for  moderate  distances  are  unified  to  one  correlation  mode  which  is 
valid  for  the  entire  range  of  correlation  distances.  With  respect  to  three- 
dimensional  cases,  introducing  remarks  are  given. in  the  appendix. 
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